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Abstract. In this manuscript, we formulate the problem of source localiza¬ 
tion based on Time Differences of Arrival (TDOAs) in the TDOA space, i.e. 
the Euclidean space spanned by TDOA measurements. More specifically, we 
show that source localization can be interpreted as a denoising problem of 
TDOA measurements. As this denoising problem is difficult to solve in gen¬ 
eral, our analysis shows that it is possible to resort to a relaxed version of it. 
The solution of the relaxed problem through linear operations in the TDOA 
space is then discussed, and its analysis leads to a parallelism with other state- 
of-the-art TDOA denoising algorithms. Additionally, we extend the proposed 
solution also to the case where only TDOAs between few pairs of microphones 
within an array have been computed. The reported denoising algorithms are 
all analytically justified, and numerically tested thorough simulative campaign. 
TDOA space and TDOA denoising and TDOA redundancy and Source local¬ 
ization 


1. Introduction 


Source localization is a research theme that has signihcantly grown in popularity 
in the past few decades, and whose interest ranges from audio to radar. As far as 
audio signal processing is concerned, several applications including teleconferencing 


(D’Arca et al, 20141, audio-surveillance (Valenzise et al 2007) and human-machine 


interaction (Trifa et al 2007) can benefit from the knowledge of the source loca¬ 


tion. Among the techniques that are available in the literature (Benesty and Huang 


2004), those based on Time Difference Of Arrival (TDOA) measurements are partic¬ 
ularly appreciated for their modest computational requirements. TDOAs, in fact, 
are usually estimated through peak-picking on the Generalized Cross Correlation of 


the signals acquired at microphone pairs (Knapp and Carter, 1976 lanniello 1982), 


E-mail addresses: {marco. compagnoni, antonlo. canclini, paolo .bestagini , 
fabio.antonacci, augusto.sarti, stefano.tubaro jOpolinii.it. 


1 




























2 


M.COMPAGNONI, A.CANCLINL P.BESTAGINI, F.ANTONAGCI, A.SARTI, S.TUBARO 


or on the whole set of microphones (Hu and Yang 2010 Chen et ah 20021. TDOAs 


can be easily converted to Range Differences (RD), once the sound speed is known. 
The source location is then found as the point in space that best fits the RD mea- 


surements according to properly defined cost functions ( 

Hahn and Tretter 

1973 

Stoica and Nehorail|19881|Schau and Robinson, 19871 Huang et all 

2001 

Bee 

c et al 


2008 Schmidt 19721. More recently, the widespread diffusion of sensor networks 
stemmed an interest in source localization also in other research communities, such 


as remote sensing and radar ( 

Koch and Westphal 

1995 

Yimin et al 

2006 

). In this 

context range differences are 

obtained from TDOAs ( 

<)ehu et al 

2009) 

, or from 


energy measurements (Hu and Li, 2002). 


The main drawback of TDOA-based localization techniques lies in their sensi¬ 
tivity to measurement noise. In particular, we can distinguish between additive 
noise (generally due to sampling in the time domain, circuit noise, but also other 
physical phenomena) and outlier measurements (produced by reverberation or in¬ 
terfering sources). Outlier identification and removal has been widely studied in 
the literature (see for instance ( Scheuing and Yan^ [2008 ) and references therein; 
(Canclini et al, 2013) and (Canclini et al, 2015| )). Therefore, applying one of these 
techniques it is possible to remove outliers from the pool of available measurements. 
Nonetheless, additive noise still impairs the localization accuracy. 

In this manuscript we interpret the problem of source localization studying the 
effect of additive noise on TDOA measurements using the TDOA space formalism, 
i.e., a space in which a set of measured TDOAs is mapped into a point. The sensi¬ 
tivity to noise afflicts also Range Differences obtained from energy measurements. 
Indeed, the hostile propagation conditions yield a difference of the measured energy 
from the ideal free-field assumption. In the following we will specifically refer to the 
problem of localizing acoustic sources, but the theory can be readily applied also 
to other kinds of signals. The concept of TDOA space is not novel and was first 
introduced in ( [Spencer 2007) for localization purposes. From that representation, 
a TDOA map (from the space of source locations to the space of TDOAs) was later 
introduced and analytically derived in (Compagnoni et al 2014), which proposed an 


exhaustive analytic study of the identifiability and invertibility of this map for the 
three-microphone case. In the most general case, given a set of n -|- I microphones, 
q = n{n -|- l)/2 TDOAs can be computed considering all the possible microphone 
pairs. In a noiseless scenario, however, we can always find an independent set of n 
such TDOAs that we can compute all the other TDOAs from. This is why most 
TDOA-based algorithms define a reference microphone, with respect to which the 
n independent TDOAs are computed. In the TDOA space, this corresponds to the 
fact that TDOAs lie on a linear subspace Vn of the q-dimensional TDOA space. 
This subspace can be computed in closed form through simple considerations. Fea¬ 
sible TDOAs (i.e. points in the TDOA space that correspond to source locations) 
are bound to lie in a region 0„ C R„. In (Compagnoni et al 2014) authors derive 
On in terms of real algebraic geometry. 

Working in the TDOA space essentially means solving an estimation problem 
in its dual space. As typically done in estimation theory, using a dual domain 
enables to split a problem into two parts. In our case, using the TDOA space 
formalism, source localization can be interpreted as a two-step procedure: i) a 
denoising operation, which consists in removing or attenuating part of the additive 
noise; ii) the application on the denoised TDOAs of a simple mapping from the 
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TDOA space to the geometric space. Starting from this perspective, in this work 
we provide a deeper investigation on the geometrical characteristics of the TDOA 
space. More specifically, we first derive the correct denoising formulation that fully 
describes the source localization problem. As the denoising problem formulated this 
way is not easy to deal with, we resort to a relaxed version of it, which exploits the 
linear subspace 14,- In particular, we show that additive noise can be decomposed 
into the sum of two orthogonal components, and the relaxed problem formulation 
aims at reducing only one of them, still positively impacting on source localization. 


This relaxed version of the problem was implicitly solved in (So et al 20081 and 


(Schmidt 1996), where the authors derive closed-form expressions for converting 
the full set of TDOAs to the nonredundant one. However, authors in (So et al 


2008) and (Schmidt 1996) limited their analysis to simulations showing that this 


conversion is able to reduce the impact of noise in localization accuracy. Instead, 
working in the TDOA space paves the way to a deeper understanding of the impact 
of relaxed denoising on source localization. In particular we will: 

(1) analytically prove the positive effect of denoising on source localization 


through a set of solid theorems, thus also theoretically validating (So et al 


2008) and (Schmidt 1996); 


(2) find a solution to the relaxed denoising problem when some TDOA mea¬ 
surements are not available; 

(3) quantify analytically the improvement in localization accuracy brought by 
the use of the denoised TDOAs given a specific localization algorithm in 
use. We accomplish this analysis by means of the error propagation theory 


introduced in (Compagnoni et al 2012). 


We test the presented algorithm also under different noise hypotheses to show 
that it works also if the underlying assumptions are not strictly verified. In particu¬ 
lar, Monte Carlo simulations were carried out to show how different state-of-the-art 


techniques (the SRD-LS algorithm (Beck et al 2008), Least Squares (Smith and 


Abel, 1987)) and Gillette-Silverman (Gillette and Silverman[ 2008)) methods ben¬ 


efit from denoising, approaching the RMSE Lower Bound (RLB) implied by the 
Gramer-Rao Lower Bound (CRLB) (Benesty and Huang 2004). We also show 


that it is possible to perform denoising on a set of TDOAs including q — s, s > 0 
measurements, with an apparent advantage in terms of accuracy. 

The rest of the manuscript is structured as follows. In Section|^we deeply intro¬ 
duce the formalism of TDOA space. In Section]^ we interpret the problem of source 
localization within this context. In Section]^ we provide the denoising formulation 
of the source localization problem, also reporting the relaxed problem version and 
an algorithm for its solution. A parallelism with related state-of-the-art works is 
also provided. In Sectionwe analytically prove the positive impact of denoising on 
source localization, and provide additional simulative analysis. Section|^is devoted 
to denoising problem formulation and solution in case some TDOAs are missing 
within the pool of measurements (i.e., we measure TDOAs using only a few pairs 
of microphones). Finally, Section remarks some final conclusions highlighting 
possible open future research lines. 


2. Theoretical background 


In this section we offer the reader some background that will simplify the reading 
of this article. In particular, we first provide the formal definition of the TDOA 
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space. Then, we give the interpretation of noisy measurements and source localiza¬ 
tion problem in the TDOA space. 

2.1. The TDOA space. The ideas of the TDOA space; the feasible set of TDOA 
measurements; and the TDOA map appeared in several manuscripts concerning 


multilateration, see for example (Schmidt 1996 Grafarend and Shan 2002 Spencer 


2007 Compagnoni et al, 2014). These concepts are the essential ingredients for the 


mathematical definition and analysis of many problems involving TDOA measure¬ 
ments, such as source localization, synchronization and calibration of the receivers. 


A recent example in this direction can be found in (Alameda-Pineda and Horaud 


2014), where the TDOA space formalism is used for defining a novel algorithm 


to estimate the TDOAs and concurrently locate the source. In the following, we 
present the basic definitions and properties regarding the TDOA space. 

Let nii = {xi, yi, Zi)'^, i = 0,... ,n he the sensor locations and x be the source 
position in the 3D Euclidean space K^. For notational simplicity, and with no loss 
of generality, in what follows we assume the sound speed to be equal to 1, so that 
the noiseless TDOAs correspond to the range differences. This way, given any pair 
of sensors (mj, mi), n > j > i > 0, the relative TDOA is a function of the source 
position X and it can be defined as 


( 1 ) 

where 


Tji(x) 


(2) rj,(x) = ||x-mjll - ||x-mill. 

If we collect the q = range differences in a ^-dimensional vector, we obtain 

the map 


( 3 ) 


(rio(x),r 2 o(x),... ,r„„_i(x))^ 

has been called the complete TDOA map, while 
‘ is referred to as the TDOA space or r-space. 


In (Compagnoni et al 2014), r* 
the resulting target set of x, 

Clearly, a point in the TDOA space corresponds to any set of TDOA measurements. 
Moreover, in a noiseless scenario, the subset of the r-space containing the TDOAs 
generated by all the potential source positions coincides with the image Im(T^) of 
the TDOA map, and we call it 0„. This means that any collection of noiseless 
TDOAs defines a point t = (tiq, ..., r„„_i)’^ G 0„ and viceversa. 

2.2. The 2D case with n=2. The study of the properties of is a fundamental 
step towards a deeper understanding of the geometrical acoustics model for TDOA- 
based localization. However, since its inherent complexity, the full description of 
the general case of r* goes beyond the scope of this manuscript. In this section, 
we summarize the main results contained in ( [Compagnoni et~S 2014 Compagnoni 
and Notari 2014) on the minimal case of two dimensional source localization, with 
three synchronized and calibrated sensorsj^ We report this analysis because this 
is the minimal non trivial case of TDOA-based localization, the only one that has 


^In order to simplify the presentation, we consider only the case with the microphones in 
general position on the plane, i.e. they do not lie on a line. The interested reader can find the 
complete analysis for every scenario and the proofs in the original manuscripts. 
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been exhaustively studied and where one may observe some important features 
characterizing every localization model. We will return on this model in Section]^ 
In the planar case, the set 02 is a surface embedded into being the image of 
the restriction of the TDOA map to (by abuse of notation, we continue to 
name it t^)- Actually, one can interpret as a (radical) parameterization of 02- 
Moreover, it is well known that the three TDOAs are not independent, since they 
satisfy the zero-sum condition (ZSC) ( Scheuing and Yang| 20061. Indeed, the linear 
relation 7-21 (x) = 7-20 (x) — rio(x) holds for each x € Geometrically speaking, 
this means that three noiseless TDOAs are constrained on the plane 


(4) 


1^2 = {t e I Tio - 720 + 721 = 0} C ■ 


and so 02 C V 2 - 

Because of the above linear relation, in the literature it is customary to work 
with a reference microphone, for example mo, and to consider only the two TDOAs 
7io(x),72 o(x). Mathematically speaking, let us define the reduced TDOA map 

72: K2 ^ ^2 

X -^ (710 (x), 720 (x)) 


and let us consider the projection map pa : forgetting the third coordinate 

721 of the 7-space. Then, we have T 2 = p^oT^ andpa is a natural bijection between 
1111 ( 72 ) and Im(72). Hence, one can investigate the properties of the noiseless TDOA 
model by studying the simpler map 72 . For the sake of simplicity, in Figure we 
draw Im(72) for the configuration of the microphone at mo = (0,0)^, mi = 
(1, 0)^ and m 2 = (1,1)^. Symbols defined therein are introduced in the next few 
paragraphs. Figure]^ shows its relation with Im(72 ). 

Let us define the displacement vectors dji = mj — m^; their Euclidean norms 
dji = IIdjiII, i,j = 0,1,2; the scalar W = det ( dio d 2 o ); and the matrix 


H = 


0 

1 


-1 

0 


First of all, 1111 ( 72 ) is contained into the hexagon P 2 defined by the triangle inequal¬ 
ities: 

{ —dio < 7io < dio 
— d20 < 720 < ^20 
— d21 < 720 ~ Do < d2l 

In particular, the vertices = {dio,d 2 o), = (—^ 107^21 — dio): dd? = (^21 — 
<^ 207 —d 2 o) of P 2 correspond to the pairs of TDOAs associated to a source at 
mp, mi, m 2 , respective ly. Then, by following the analysis contained in Section 
6 of (Compagnoni et al 20141, for any 7 = ( 710 , 720 ) G we define the vectors 

(7) 

idlo - 7 |o) dio - (d?o - D^o) ^20 


v( 7 ) = H (720 dio - 710 d2o) 7 
and the polynomials 
( 8 ) 0(7) = ||v( 7 ) 


lo(7) =H- 


2W 




6(t) = v(7)^ • 10 ( 7 ), c{t) = ||lo(7)|l 


Im(72) and the admissible source positions 72~^(7) can be computed in terms of 
these polynomials. 

We have the following facts: 
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Figure 1. The image of T 2 is the gray subset of the hexagon P 2 with 
continuous and dashed sides. In the light gray region E~ the map 
is 1-to-l, while in the medium gray region [/o U C/i U U 2 the map T 2 is 
2-to-l. Let us observe that [/oUl/iUt /2 C C'^nP 2 . The continuous part 
of the boundary of the hexagon and the blue ellipse E, together with 
the vertices P*, are in the image, and there T 2 is 1-to-l. The points 
and the dashed boundaries do not belong to Im(T 2 ). 


• a{T) — 0 defines the unique ellipse E tangent to every facet of P 2 - We name 
E~ and E~^ the interior and the exterior regions of E where a(T) < 0 and 
o(t) > 0 , respectively; 

• 6(r) = 0 defines a cubic curve C. We name C~ and C+ the regions where 
6(r) < 0 and 6(r) > 0, respectively; 

• c(t) is a quartic non negative polynomial. 


In (Compagnoni et al 20141 it has been proved that the image of T 2 is the set 


(9) 


Im(T2) = E~ ((7+ n P 2 ) U i?°. 


For each r S Im(x 2 ) we have at most two admissible source positions, whose 
coordinates are given by the formula 


(10) x±(r) = mo + 1o(t) + A±(x)v(t), 

where A±(t) are the solutions of the quadratic equation a(x)A^+26(T)A+c(T) = 0 : 

^ -&('r) ± - a(r)c(T) 

Ai(r) =--. 

For T G {E~ U E U dP 2 ) H Im(x 2 ) we have to take only the x+(t) solution and so 
we have uniqueness of localization. On the complementary set, that is the union 
of the three disjoint sets [/qi Ui and U 2 depicted in medium gray in Figure]^ the 
map T 2 is 2-to-l and there is an intrinsic ambiguity in the source position between 
the two solutions x±(t). 
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Figure 2. The image of t* is the green subset of the hexagon 1^2 C 
V 2 , while the image of T 2 is the red subset of P 2 . There is a 1-to-l 
correspondence between Im(T 2 ) and Im(T 2 ) via the projection map pa. 
In the lightly shaded regions, the TDOA maps are 1-to-l, while in the 
more darkly shaded regions the maps are 2-to-l. 


For the sake of completeness, in Figurej^we depict the corresponding localization 
regions in the T-plane. Roughly speaking, we have the preimage of the interior 
of the ellipse E~ = where the TDOA map is 1-to-l and the source 

localization is possible, and the preimages Ui = T 2 ~^{Ui), for i = 0,1, 2, where the 
map is 2-to-l and there is no way to uniquely locate the source. The region of 
transition is the bifurcation curve E = T 2 ~^{E), that is a quintic algebraic curve 

2014[ ) consisting of three disjoint and unbounded arcs, 
one for each arc of E contained in Im(T 2 ). As a point r in one of the Ui gets close 
to A, the solution x+(x) gets close to a point on E, while x„(x) goes to infinity. 
The sets E~,Uo,IJi,U 2 are open subsets of the T-plane, separated by the three 
arcs of E. 

Finally, the union D of the six dashed half-lines outgoing from the receivers is 
the degeneracy locus of the TDOA map, where the rank of the Jacobian matrix of 
X 2 drops. D is the preimage of the six segments in dP 2 H Im(x 2 ). On D the two 
solutions x±(x) are coincident, thus the TDOA map is 1-to-l. Furthermore, D 
divide each Ui into two connected components and T 2 is a bijection between each 
of them and the corresponding Ui. 

2.3. The general case. As we said above, some of the properties we described in 
the minimal planar case are common to every TDOA-based localization model. In 
particular, we have the following proposition. 

Proposition 2.1. Let us take n + 1 sensors at mo, ■ • •, nin in where n > 2. 
Then, 0„ is a subset of the n-dimensional linear subspace C defined by 


(Compagnoni and Notari 
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Figure 3. The different localization regions and the curve E in the x- 
plane. The microphones are the marked points mo = (0,0), mi = (1, 0) 
and m 2 = (1, !)• The curve E separates the light gray region E~, where 
the map T 2 is 1-1 and it is possible to locate the source, and the medium 
gray region Uo VJ Ih VJ C/ 2 , where T 2 is 2-1 and the localization is not 
unique. On the dashed lines the localization is possible but very sensitive 
to the measurement noise. 


equations 

(11) - Tio + Tjo -Tji = 0, 0 < i < j < n , 

representing the ZSCs for all the microphone triplets containing mo. 


Proof: in a configuration with n + 1 microphones, the maximum number of in¬ 
dependent TDOAs is equal to n. In particular, if we take mo as the reference 
microphone, the n TDOAs {tio(x), ... ,rno(x)} are independent, while the others 
satisfy equations 0’ as can be easily verified using definition ([^ . These are q — n 
independent homogeneous linear equations, therefore they define an n-dimensional 
linear subspace I 4 of the r-space and 0„ is a subset of □ 


This means that, also in the general case, the set of feasible TDOAs is contained 
into an n-dimensional linear subspace of the TDOA space K®. This property stays 
at the basis of the denoising procedure that we describe in the next sections. How¬ 
ever, 0„ is strictly contained in I 4 . Indeed, since 0„ is the image of through 
the almost everywhere smooth function t„, its dimension is equal to 3. As above, 
one can consider Tn as a radical parameterization of 0„. This way the feasible 
set becomes a topological manifold (possibly with a boundary, as for 02 ), that is 
almost everywhere differentiable. Moreover, we can reasonably conjecture that 0„ 
can be described again in terms of algebraic equations and inequalities, hence it is 
a so called semialgebraic variety (iBasu et all 12006 1. 


3. Interpretation of source localization in the TDOA space 


The TDOA space formalism that we just introduced can be used to provide a 
geometric interpretation of TDOA-based source localization problem. As a mat¬ 
ter of fact, other geometric interpretations have been proposed in the literature 
(Bestagini et al 20131. However, working in the TDOA space also highlights that 
source localization can be solved as a TDOA denoising problem. In this section. 
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we provide such interpretation of source localization in the TDOA space, starting 
from a commonly used statistical noise model. 

3.1. Statistical noise model. In the presence of measurement errors, we must 
resort to statistical modeling. In this manuscript we assume the TDOAs associated 
to a source in x to be described by 

(12) f*(x) = T*(x)+e, 

where e ^ N{0,'E) is an additive Gaussian noise. Also techniques that compute 
Range Differences from other measurements, such as energy, are prone to additive 
noise. In this latter case, in particular, the magnitude of the additive noise becomes 
relevant. Under the assumption in (|12p, the probability density function (p.d.f.) of 


the TDOA set is (Benesty and Huangl 20041 


dill), 


(13) 


p('r;<(x),s) = 




e 2 


v/(2^)«|S| 


Since the covariance matrix S is symmetric and positive d ehned, from a geometric 
standpoint (see, for example, (Amari and Nagaoka 20001) the Fisher matrix 
defines a scalar product on 


(14) 


(Vl,V2)s~i = ^2, Vi,V2e 


This way, the TDOA space turns out to be equipped with a Euclidean structure, 
whose distance is known in the statistical literature as the Mahalanobis distance: 


(15) 


.-1 = 


V e 


With this setting, the p.d.f. (13) can be rewritten as 

-UU-<(x) 


(16) 


p('r;<(^)>^) = 


V(2^)«|S| ’ 

which depends only on the Mahalanobis distance between r and (x). 

3.2. Source localization. The first application of the TDOA space and map was 
in the study of the TDOA-based source localization (see (Spencer 2007 Com- 
pagnoni et al 20141). As a matter of fact, the fundamental questions in localiza¬ 


tion problems can be readily formulated in terms of t*. In a noiseless scenario, the 
analysis of the existence and uniqueness of localization is equivalent to the study 
of the set 0„ and the invertibility of As we saw in Section 2.2 for the minimal 
planar case, for a given r G there exists a unique source at position if, 

and only if, x is a point lying on a region of 0„ where the TDOA map is 1-to-l. 

In case of noisy measurements, the data errors force the point f not to lie on 0„. 
Therefore, to localize the source one needs an estimation procedure. Most of the 
algorithms proposed in the literature rely on the optimization of a cost function, 
based on some criterion that can be either statistically motivated (e.g.. Maximum- 


Likelihood estimation ( 

Benesty and Huang 

2004 

l)or 

not (e.g., linear Least Squares 

(Smith and Abel 

1987 

Gillette and Silverman |2008 

I, Squared Range-Differences- 

based least squares estimation ( 

Beck et al 

20081, etc.). The source position is thus 


found as the point x that minimizes a suitable non-negative cost function /(t,x), 
defined so that its value is zero in noiseless conditions, i.e.. 


( 17 ) 


./«(x),x) = 0 
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for every x S In mathematical terms, source localization is therefore formulated 
as 


(18) 


X = argmin /(t, x) 

xeR3 


In the TDOA space, the estimated source position is associated to the feasible 
point f = S 0„. Therefore, any localization algorithm maps a noisy TDOA 

vector T onto a feasible TDOA vector r. It is worth noticing that different al¬ 
gorithms may produce different source estimates, corresponding to mappings onto 
likewise different feasible TDOA vectors. A special case occurs when the input 
TDOA vector is feasible, i.e., when r S On- In this case, there exists a point x so 
that T = T^(x). Thus, in virtue of we expect any algorithm to produce the 
same estimate x = 

Let us now focus on the Maximum-Likelihood (ML) estimator, which is known to 
be optimal in the statistical sense. For the Gaussian noise model described above, 
the ML localization problem can be formulated as ( [Benesty and Huang 2004) 

(19) 


Xml = argmax 

xG 


By defining /ml(t, x) = ||t - t*(x) 
( 20 ) 


R3 (27r)9|S| 

we have 


Xml = argmin /ml(t,x). 
xeR3 


In the TDOA space framework, the ML estimator has a neat geometric inter¬ 
pretation. Indeed, solving the ML problem is equivalent to finding the point 
"Tml = t^(xml) G On at minimum Mahalanobis distance from r: 

(21) Tml = argmin IIT - 

TGSn 

It trivially follows that the source position x estimated by means of a generic 
localization algorithm leads to a point T = (x) € 0„ such that 

(22) ||t-t|1s-i > ||t-Tml|Is-i ■ 


In other words, the distance HF —T|j^.-i is bounded from below by the one obtained 
through ML estimation. 

In the light of the above considerations, we can interpret source localization in 
the TDOA space as a two-step procedure: 

(1) mapping of the noisy TDOAs t onto the corresponding feasible vector 
r € 0„, accordingly with the chosen optimization criterion; 

(2) recovering of the estimated source position as x = 'r^~^(T). 

Note that step 1 represents a TDOA denoising operation. Moreover, once this step 
has been accomplished, step 2 is straightforward when the mapping r* is 1-to-l, 
that is the standard case for a sufficiently large number n of sensors in general 
positions. From this perspective, source localization can be considered as a TDOA 
denoising problem. 


4. Denoising of TDOAs 

As discussed in Section source localization can be solved in the TDOA space 
as a denoising problem. Given a noisy TDOA vector t, this corresponds to finding 
a feasible (i.e., denoised) TDOA vector belonging to 0„ according to some criteria 
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that depends on the chosen cost function /(t,x). If we consider the Maximum- 
Likelihood formulation of the problem, the criterion is readily formulated as in 
(211. Therefore, statistically speaking, the best achievable denoising corresponds 


to finding the feasible TDOA vector at the minimum Mahalanobis distance from 
the noisy one. 

Despite the error function ||t — could suggest a standard weighted linear 


least-squares solution (Teunissen 2000), the problem in (21) can not be solved in 


a linear fashion. Indeed, the search space is 0„, which turns (21) in a difhcult, in 


general non convex, optimization problem. However, in the following we investigate 
the possibility of relaxing the problem in (21), leveraging on the fact that 0„ is 


contained in the linear subspace 14, of the TDOA space. 


4.1. Prom the complete to the relaixed denoising problem. We claimed 
above that solving the ideal denoising problem (21) is a non trivial task. How¬ 


ever, by working in the TDOA space we can subdivide such problem in two distinct 
steps. Let us consider the orthogonal projection 7 ^(t;S) of t e on 14, with 
respect to the scalar product ( , )s-i. From the general properties of this map, we 
have the following equivalences: 


Tml — 


s- 


argmin ||t — t\ 

argmin|l(f-lP(T; S)) (P(t; i:)-r)|||. 

-rG0„ 

argmin(||f-P(f;S)|| 


2 


|lP(f;S)-i 




= argmin ||7^(r;S)-T| 




S-1 ’ 


where in the third equality we used the fact that r — P(t; S) and 'P(t; S) — t are 
orthogonal each other, thus 

(f - 'P(t; S), 'P(f; S) - Tml)s-i = 0 . 

This means that the ML estimation gives exactly the same results if we start from 
the original data f or the projected one P(t; S). 

We explicitly show this fact in Figure]^ for the case of planar localization with 
n = 2. In this case one has two different situations, exemplified for the measure¬ 
ment points Ti and T 2 , respectively. In the first case, the orthogonal projection 
P(fi;S) e 02, thus it coincides with its ML estimate Tml,i- Differently, for the 
second point one has P(t 2 ;S) ^02- In this case, finding the corresponding ML 
solution Tml .2 corresponds to finding the closest point to P(t 2 ;S) in 02 , which 
implies solving a complicated minimization problem in 14. Indeed, as we explained 
in Section |2.2[ the feasible set 02 has a non trivial structure. In particular, its 
boundary is the union of six segments and three arcs of ellipse. If we define Tml,2 
as the closest point to T 2 lying on one of these sets, then it is necessary to develop 
an ad hoc algorithm for finding it. Another complication is that 02 is not a closed 
set. This implies that if rML ,2 lies on the ellipse, then it does not correspond to 
a true source position because it is not part of 02. Conversely we should consider 
'Eml ,2 as the TDOAs associated to a source placed at infinity. 

Thanks to the results and discussion contained in Section [2.3[ we can generalize 
the previous analysis. The denoising problem (21) can be subdivided into two 
subproblems. 

















12 M.COMPAGNONI, A.CANCLINL P.BESTAGINI, F.ANTONAGCI, A.SARTI, S.TUBARO 



Figure 4. Denoising in the TDOA space, for the case of planar local¬ 
ization with n = 2. A generic set of noisy TDOAs t do not lie on 02. 
The ML estimation finds the closest point tml G 02 to t. The solution 
of the relaxed denoisig problem is 'P(t, S) G V 2 . For the point ti it 
lies on 02, then the projection V{f\, S) coincides with the ML solution 
Ti. For the point T 2 this is not true, thus P(ti, S) ^ f 2 - In any case, 
the projection gives a better estimate of the noiseless measurements t, 
closest to Tml with respect to the original data point r. 


( 1 ) 


( 23 ) 


The easiest part is the projection on the linear subspace Vn 
the relaxed denoising problem 


We can call it 


7^(t; S) = argmin 

TG Vn 


T — T 




in comparison to the complete denoising problem (211. Being the search 


set in (23) a linear subspace, the problem admits a closed solution. 


(2) The hardest part is the projection of V{f\ S) onto 0„ 


Tml = argmin ||'P(t; S)-t|||,_i , 

From the discussion contained in Section |2.3[ the difficulties are twofold: 
we have not the analytic description of 0„ and, in any case, the feasible set 
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is a complicated three dimensional semialgebraic variety embedded in the 
n- dimensional linear subspace Vn- 

The previous observations are the original starting points for our interpretation 
and for the statistical justification of the relaxed denoising procedure. Indeed, 
observation 2 confirms the well known fact that solving ML estimation is a too 
difficult task, due to the non-linearity and non-convexity of the feasible set 0„. 
However, observation 1 states that ML is composed by two parts with different 
difficulties. We can easily envision that the projection of the measured TDOAs onto 
the linear subspace Vn leads to improvements in terms of localization accuracy, and 
this will be confirmed analytically in the forthcoming sections. 

4.2. The relaxed denoising algorithm and its statistical analysis. We can 

now proceed with a rigorous formulation of the intuition discussed above and with 
a precise definition of the relaxed denoising algorithm. 


4.2.1. The projection as a sufficient statistic. 

Theorem 4.1. The orthogonal projection o/t G K® on Vn, with respect 

to the scalar product { , a sufficient statistic for the underlying parameter 

X. 


Proof: Let us start from 




= \\f-V{f-,'E)\ 


||7^(f;S)-<(x) 


1 ^ 

s- 


Therefore, the probability density function (13) can be rewritten as 








Then, the proof follows as a consequence of the Fisher-Neyman factorization theo¬ 
rem (Lehmann and Casella 19981. □ 


Theorem |4.1| states that the component of t orthogonal to the linear subspace 14, 
does not carry information on the source location and that we can remove it without 
corrupting the data. 


4.2.2. The algorithm. For any r G M'^, the projection 'P(t; S) can be computed in 
closed form in two steps: 

1. If we group all the equations 0’ we end up with the homogeneous equation 
system 


(24) Cx = 0, 

where C is the {q — n, q) matrix 


( 

1 

0 

0 0,-1 0 •• 

■ 0 \ 

-1 

0 

1 • • 

0 0 1 0 -1 •• 

. 1 . 

0 

\ 0 

0 

0 •• 

• -1 1 ' 0 0 •• 

• -1 / 


The solution of this linear system is 14 = ker(C). In particular, we can find an 
orthonormal basis {vi,..., of 14, if necessary by using Gram-Schmidt 

algorithm defined with respect to the scalar product (14). 
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11. The projection map V is defined on f € as 

(25) = (f,vi)s-i vi + ... + (f,v„)s-i v„. 

Let Cji, n > _) > i > 1 be the vectors in the standard basis Bq of With 
respect to Bq, the projection is represented by the {q,q) matrix 

(26) P=[iP(eio;S) ... S) ] . 

Consequently, the set of denoised TDOAs is obtained as 

(27) iP(f;S)=PT. 


4.2.3. The analysis of the noise reduction. We can now compute the noise reduction 
on the TDOAs due to the relaxed denoising procedure defined above. Preliminarily, 
we prove the following Lemma. 


Lemma 4.2. The covariance matrix S defines a Euclidean structure on and 
the matrix P^ represents an orthogonal projection with respect to the scalar product 
( ) )s- 


Proof: By construction S is symmetric and positive defined, therefore it defines a 
scalar product on K'?. To prove that P^ represents an orthogonal projection with 
respect to this Euclidean structure, we have to show that P^P^ = P^ and SP^ = 
PS. On the other hand, we know that the matrix P represents an orthogonal 
projection with respect to ( , )s-i, hence PP = P and S^^P = P^S~^. By using 
these identities, we have: 

. P^P^ = (PP)^ = P^, 

• PS = ss-ips = SP^S'^S = SP^, 


which completes the proof. 

A consequence of Lemma 4.2 is that 


!i: > 


for any v G 


□ 

This is 


useful for proving the next Theorem. 


Theorem 4.3. Let S be the covariance matrix of t. Then: 

(1) the covariance matrix o/P(t;S) is S' = PSP^; 

(2) we have S ^ S', i.e. S — S' is positive semidefinite. 

Proof: the first claim follows from the general transformation rule of the covariance 
matrix under linear mapping of t. Using this fact and Lemma |4.2[ we have that 
v^(S — S')v = |jv|||. — ||P^v|j|. > 0 for any v S M'?, which completes the proof. 
□ 

Theorem |4.3| states that the relaxed denoising procedure always reduces the noise 
on the TDOA dataset and it gives the way to quantify such reduction. 


4.3. Relation to state-of-the-art algorithms. As already mentioned in the In¬ 
troduction, the algorithm presented to solve the relaxed denoising problem is a dif¬ 
ferent interpretation of other methods proposed for reducing the noise on TDOAs, 


exploiting data redundancy. In particular, authors in (Schmidt 1996) and (So et al 


20081 use the constraints in (111 to relate the full set of q TDOAs r to n nonredun- 
dant TDOAs referred to a sensor. Selecting the first sensor as the reference one. 
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the nonredundant set is = (tiq , T20 , ■ • ■ , t„o) 
the two sets is given by 


The linear relation between 



fl 1 


/ -1 1 0 

-1 0 1 

0 0 \ 

0 0 

T = Gtnr , with G = 

Y 

where Y = 

\ 0 0 0 •• 

.-11/ 


and I„ is the identity matrix of order n. Given the measured TDOAs t, in (Schmidt 


19961 the nonredundant set is estimated in the least squares sense as 

= (G^G)-iG^f . 


"^NR 


In (So et al 20081 this result is generalized accounting for the covariance structure of 


noise. To this purpose, the noise covariance matrix S is introduced in the weighted 
least squares solution 

(28) Tnr = (G^S~^G)-iG'^S~^f . 

Note that the two procedures coincide if S = cr^I. 

The nonredundant TDOAs t^r can be considered as the denoised measurements 


corresponding to the original set f. As a matter of fact, Tnr, computed as in (28), 


coincides with the first n components of the orthogonal projection V{t; S). At this 
respect, we can be more precise. 

Proposition 4.4. Given the Euclidean structure on defined hy the matrix 
P that represents the projection map V, with respect to the canonical basis Bq, is 

P = G(G^S~^G)-iG^S'^ 

It follows that Tffji is a sufficient statistics for f. 

Proof: First of all, we show that Im(P) = I4ij^ We begin by observing that the 
matrix G^S“^G is invertible. Indeed, its kernel is given by the vectors x G K" 
such that G^S~^Gx = 0. In this case, we have 


'^G^S-iGx= ||Gx|l|.i =0, 


hence Gx = 0. Being rank(G) = n, it follows that ker(G^S ^G) = ker(G) = {0}. 
From the previous point, we have 

rank((G^S-iG)-iG'^S-i) = rank(G^) = n 

and so Im((G^S-^G)-^G^S-^) = K". It follows that Im(P) = Im(G) = I 4 . 

Now, it is straightforward to check that P^ = P and S~^P = thus P 

represents the projection map on Im(P) = Vn- 

The last statement follows from the I-to-I relation V(t; S) = Gtnr- □ 

In the light of the above considerations, all our theorems and simulative results 


can be considered as further validations of state-of-the-art works in (Schmidt 19961 


and (So et al 2008). In the following sections we will prove, both analytically and 


numerically, that any source localization algorithm benefits from denoised TDOAs. 
Moreover, we will see that this holds also when the full TDOA set is not completely 
available. 


^With slight abuse of notation, in the proof we identify a matrix with the associated linear 
map defined through matrix multiplication. 
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5. Impact on source localization 

To this point we have shown how it is possible to reduce the noise on a TDOA set 
T, thus approaching Tml- The goal of this section is to prove that source location 
estimate benefits from the use of a denoised TDOA set. To this purpose, we first 
provide a solid theoretical analysis that demonstrates our claim. Then, we further 
investigate the effect on source localization by means of an extensive simulative 
campaign. 


5.1. Impact on sub-obtpimal localization algorithms. The present analysis 
is mainly based on the results in Section [4. 2. 3[ In great generality, let us consider 


a given localization algorithm based on the minimization of a certain cost function 
/(x, c), where c are the input TDOA data. Theorem 4.3 allows us to compare 


the accuracy of x = arg min/(x,T) and x' = arg min/(x, 7 ^(t; S)). Indeed, the 
first order error propagation formula given in ( Compagnoni et al[ 20121 relates the 
covariance matrices S, S' to Sx, S^, respectively: 


(29) 


Sx = A(x)SA(x)^ and S(j = A(x)S'A(x)^ 


where A(x) is defined by equation (26) from (Compagnoni et al 2012). Then, by 


easily adapting the proof of Theorem |4.3[ we end up with 
Corollary 5.1. at first order approximation, Sx ^ S(j. 


We remark that Corollary 5.1 is valid for every choice of the cost function /(x, c). 


In particular, it includes the special cases of /(x, c) explicitly depending only on n 
TDOAs, the ones relative to a reference microphone. In the following paragraph 
we will report concrete examples of these facts. 


5.2. Numerical examples. We now analyze the effect of relaxed denoisi nj^ 
simulated data. Some results were reported in (Schmidt 


on 


Schmidt 

1996 

So et al 

2008 


In particular, in (So et al 2008) authors proved that the denoising procedure can 


be used to proficiently exploit the redundancy of the full set of TDOAs measured 
among all microphone pairs. Indeed, they numerically verified that ML source 
localization using the denoised TDOA set (denoted as optimum nonredundant set 
in (So et al, 2008)) meets the Cramer Rao Lower Bound (CRLB) ([Benesty and 


Huang 2004) computed considering the full set of measurements. Note that the 


analytical proof of this fact is a direct consequence of Theorem |5.1| 

Here we aim at providing a more comprehensive analysis, testing the effect of 
denoising on several source localization algorithms. In particular, we are interested 
in sub-optimal algorithms admitting an exact solution, i.e., for which the global 
optimum can be computed with no approximations. This feature is desirable in 
practice, as it prevents the risk of getting trapped in local minima, at the expense 
of obtaining a solution that does not meet the CRLB. To this end, we consider the 
following algorithms: 


• LS: unconstrained linear least squares estimator (Smith and Abel 1987), 
admitting a closed-form solution. It uses the nonredundant TDOA set, 
measured considering a reference sensor; 


^For the sake of compactness, throughout this paragraph we will use the word denoising 
referring to relaxed denoising. 
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acoustic 

source 


Figure 5. Simulation setup. 


SRD-LS: squared-range-difference least squares estimator (Beck et al 20081, 
a constrained version of the LS one, whose exact solution is computable ef¬ 
ficiently. It is based on the nonredundant TDOA set; 

GS: Gillette-Silverman algorithm (Gillette and Silverman 20081, an ex¬ 


tension of the LS algorithm accommodating the case of multiple reference 
sensors. It can be applied to the full TDOA set. 


Denoising was tested considering TDOAs corrupted by both Gaussian and non- 
Gaussian noise. The first case constitutes an ideal scenario, in which the noise 
model exactly meets the theoretical assumptions used for deriving the denoising 
procedure. The second case represents a non-ideal condition, useful to assess to 
what extent denoising is applicable when the assumptions are not obeyed. 


5.2.1. TDOAs corrupted by Gaussian noise. We simulated a compact cross ar¬ 
ray composed by 7 microphones in positions (0,0,0)^, (±0.5, 0,0)^, (0, ±0.5,0)^, 
(0, 0, ±0.5)^ m. More than 500 sources are homogeneously distributed on a sphere 
centered at (0,0, 0)”^, whose radius d ranges from 0.5 m to 2.5 m. The simulation 
setup is sketched in Figure For each source position, we computed the full set of 
<7 = 21 theoretical TDOAs t. We corrupted the vectors t with / = 5000 realiza¬ 
tions of i.i.d. zero-mean Gaussian noise with standard deviation cr, leading to the 
noisy TDOAs t^. The covariance matrix thus resulted in S = cr^Iq, where Ig is the 
identity matrix of order q. Monte-Carlo simulations were carried out considering 
the range a G [0.5 cm, 5 cm]. The corresponding denoised TDOAs 'P(Ti;S) were 
computed using (27). 

As a preliminary test, we computed the mean and the standard deviation erg 
of the residual error €i — Ti — t left on denoised TDOAs. We first observed that 
Ttg is always negligible compared to erg, meaning that the denoising procedure does 
not introduce relevant bias on TDOAs. Figure shows the standard deviation 
of denoised TDOAs erg as a function of cr, averaged for sources located at a fixed 
distance d = 1.5 m. Similarly, Figure [6b| shows erg for different distance values, fixing 
cr = 1.5 cm. In both the cases, we notice that erg (solid line) is always significantly 
below the value of the standard deviation of the injected noise cr (dashed line), 
thus confirming the effectiveness of denoising. We also observe that erg « Icr, 
independently from the source position. 
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Figure 6. TDOA residual error before and after denoising, as a function 
of: injected noise a (a); source distance d (b). 


The localization performance was then evaluated in terms of Root Mean Square 
Error (RMSE), computed as 


(30) 


RMSE(x) = 


\ 


jtl 

i=l 


where x is the nominal source position and is its estimate at the fth Monte-Carlo 
run. The three algorithms were fed with both the measured TDOAs Ti and the 
denoised TDOAs 7^(fi;S). In particular, for LS and SRD-LS we considered only 
the n = 6 TDOAs measured with respect to the first sensor, selected as the reference 
one. For GS, we considered the full set of q TDOAs, before and after denoising. 
Results are reported in Figure ^ As done before, we report results as a function 
of a when d = 1.5 m (Figure 7a); and varying the distance d when a = 1.5cm 
(Figure[7b|). For each tested algorithm, the figures show the average RMSE achieved 
before and after the denoising of the measured TDOAs. For the sake of comparison, 
we also report the RMSE Lower Bound (RLB) implied by the CRLB. As expected, 
all the algorithms exhibit improved localization accuracy when denoised TDOAs 
are used. It is interesting to observe the different behavior of the algorithms using 
only the nonredundant TDOAs (LS and SRD-LS) and that using the full TDOA 
set (GS). Indeed, denoising produces a moderate effect on LS and SRD-LS, while 
the benefit is impressive for GS. Despite using the full TDOA set, before denoising 
GS exhibits very poor performances compared to LS and SRD-LS, especially for 
distant sources and at high amount of injected noise on TDOAs. This behavior 
is not completely unexpected, as GS extends the LS approach, which is known to 


suffer from non negligible bias (Benesty and Huang 2004). We observed that, using 


the full TDOA set, the bias increases significantly, even if the variance of estimation 
is reduced due to the availability of more measurements. We preliminarily noticed 
that the bias is greatly reduced when GS is fed with the denoised TDOA set. 
Roughly speaking, this means that denoising allows GS to effectively exploit the 
data redundancy. This positively impacts on source localization with GS after 
denoising, producing an RMSE that approaches the RLB. On the other hand, the 
redundancy is only partially exploited by LS and SRD-LS, as they rely on the 
nonredundant denoised TDOA set. 
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Figure 7. Localization accuracy before and after denoising, as a func¬ 
tion of: injected noise cr(a); source distance d (b). 


5.2.2. TDOAs corrupted with non-Gaussian noise. Considering the simulation setup 
described in the previous paragraph, we repeated all the tests injecting different 
types of non-Gaussian noise in the ideal TDOAs. In particular, we focus on the 
following noise models: 


i.i.d. zero-mean uniform noise in the range 


■ ’ 2 /, 


, which mimics the 


sampling error caused by estimating the TDOA from the a Generalized 
Cross-Correlation (GCC) function (Knapp and Carter 1976) sampled at 
fs = 8 kHz; c = 343 m/s is the speed of sound; 

a mixture of i.i.d zero-mean uniform and Gaussian noise, the former dis¬ 
tributed in the range ^ j the latter with standard deviation a = 

1.5 cm. This mixture of models is suggested in (Benesty and Huang 2004); 
i.i.d. zero-mean Laplacian noise with standard deviation a = 1.5 cm. 


The localization accuracy of the three algorithms exhibits the same trend obtained 
for the Gaussian noise model. For reasons of space, here we limit to report the 
average RMSE as a function of the source distance d. Figures [8a| [8b| and |8c| show 
the results for the three considered noise distributions, respectively. Even if the 
noise models differ from the Gaussian assumption underlying the denoising theory, 
all the algorithms improve their accuracy using denoised TDOAs. As before, this 
is especially true for GS, which exhibit the most noticeable improvement. This 
suggests that denoising can be applied when the TDOA error does not strictly 
meet the Gaussian assumption. 


6. Dealing with an incomplete set of TDOAs 

There exist many scenarios where not the whole set of TDOAs is available. For 
example, when computational cost is an issue, the computation of all the possible 
TDOAs is not feasible. In the following, we adapt our previous analysis in order to 
handle relaxed denoising also in such situations. 

Let us assume that the TDOAs ,..., }, 0 < s < g, are not available and 

let S = ..., (js,*s)} be the corresponding set of indices. In this setting, 
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Figure 8. Localization accuracy before and after denoising, as a func¬ 
tion of distance d, for different TDOA noise models: uniformly dis¬ 
tributed noise (a), mixture of uniform and Gaussian noise (b), Laplacian 
noise (c). 


the proper TDOA map is 


(31) 


‘n,S 


M3 

X 



where ^^ _ 

= (Do(x), . . (x), . . Tnn-l{x)V 

and Tj^i^ (x) means that the item is missing. As before, we define the TDOA 
space as the target set M‘^“® of r* g and the image Ini(r* g) as 0„,s. Clearly, the 
TDOA map T* g is strictly related to r*. Indeed, let us consider the projection 
PS : M"^ —)■ M'^“® that takes care of forgetting the s coordinates corresponding to the 
indices in S. Then, one has s ~ PS ° '''n ©n.S = Ps(0n), where the symbol 
o denotes the function composition operator. 

In the presence of measurement errors, we assume that the TDOAs are described 
by the statistical model 


(32) + ^s, where £s -- N{0, Eg). 

As above, the Fisher matrix defines a Euclidean structure on the TDOA 

space M'^“® and the same discussion made in Section holds in this situation. 

The definition and the analysis of the relaxed denoising procedure are very similar 
to the ones made in Section [T^ First of all, we adapt Theorem |4.1[ 

Theorem 6.1. Let us take n + \ microphones at mo, ..., mn in M^, where n > 2. 
Then Qn,s is a subset of the linear subspace Vs = ps{Vn) C M'?“®. Moreover, the 
orthogonal projection 'Ps(ts;Ss) of ts S M'J“® on Vs, with respect to { , is 

a sufficient statistic for the underlying parameter x. 


Proof: PS is a linear map, therefore Vs is a linear subspace of M"? The other 
claims follow in the same way of Proposition |2.1| and Theorem |4.1| □ 


From Theorem 6.1 it follows that dim(V 5 ) < dim(P„) = n, where dim denotes 
the dimension of a vector space. It is not difficult to prove that the equality holds 
if, and only if, the set of available TDOAs contains n independent measures, for 
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example the n TDOAs calculated with respect to a reference microphone. In this 
case, the map ps is a bijection between V and ^ 5 . This means that it is possible to 
obtain the full set of q denoised TDOAs as pg^{Vsi'rs','^s)) G V- Concretely, one 
has to use the linear equations 

(33) Tik '^jk D'l — 0, i ^ j ^ hj 

to calculate the missing TDOAs with indices in S. It is important to remark that 
this operation does not increase the noise on the dataset, as indeed it would happen 
if we apply the same procedure on the original data by calculating pg^{Ts). 


6.1. Denoising algorithm. In order to explicitly construct the projection map 
Vs, let us take the {q — s, q) matrix Is defined by removing the s rows corresponding 
to the indexes in S from the ( 9 , q) identity matrix. It can be easily proven that 
Is represents the map ps with respect to the standard basis Bq and Bq-s of 
and respectively. Then, given a generic basis {vi,..., Vn} of ker(C), from 


Theorem 6.1 follows that the set {Is Vi,..., Is Vn} spans Vs- After reducing it to an 


independent set of vectors and subsequently applying the Gram-Schmidt algorithm, 
we can finally find an orthonormal basis of Vs with respect to the scalar product 


From this point on, one can proceed exactly as done in the previous sections. 


In particular, the map Vs is defined by the analogous of the formula (25) and it is 
represented, with respect to Bq-s, by a {q — s, q — s) matrix Ps. Then, the denoised 
TDOAs are 


(34) 


Ps(fs;Ss) = PsTs • 


6.2. Impact on source localization. We summarize the main facts on the de¬ 
noising procedure in the following theorem. 


Theorem 6.2. let Sg be the covariance matrix of Ts and /(x, c) be any given cost 
function, where c are the input TDOA data. Then: 

(1) the covariance matrix of Vs{ts','^s) = PsSgPs^; 

( 2 ) Ss ^ Ss'; 

(3) at first order approximation, the covariance matrices Ss,x and Ss,x of 
the estimators x = arg min/(x, Ts) and x' = arg min/(x, P5(rs; Sg)), 
respectively, satisfy Ss,x ^ ^Sx- 


Proof: The proof is similar to the ones of Theorem 4.3 and Corollary |5.1 


□ 


6.3. Numerical examples. In this paragraph we show some numerical examples, 
devoted to investigate the effect of relaxed denoising when the full TDOA set is not 
entirely available. To this end, we refer again to the simulation setup described in 
Section |5.2.1 However, we now consider the availability of the n TDOAs referred 
to the first microphone, along with z < q — n additional TDOAs. Noisy TDOAs 
were obtained corrupting their nominal values with i.i.d. zero-mean Gaussian noise 
with standard deviation a = 1.5 cm. In this scenario, denoising was accomplished 
using (341. In particular, we built the vector ts including the n + z available 
TDOAs, and we computed the projection matrix Ps accordingly. For this test 
we considered sources at a fixed distance d = 1.5 m. We tested the denoising 
procedure considering values of z in the range from 1 and q — n — 1 = 14. For all 
the I Monte-Carlo runs, we generated all the possible combinations of z TDOAs 
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Number of additional TDOAs 


(A) 


(b) 


Figure 9. TDOA residual error before and after denoising (a) and local¬ 
ization accuracy (b). Both plots are function of the number of additional 


TDOAs 


extracted from the last n — q entries of the vector ti. As before, we considered 
the LS, SRD-LS and GS algorithms for source localization. The results, averaged 
among all the noise realizations and all the combinations, are reported in Figure]^ 
In particular. Figure shows the residual error on TDOAs after denoising, while 
Figure highlights the impact of denoising on localization. Note that we added 
to the graphs the points at z = 0 (i.e., when only the n TDOAs referred to the first 
microphone are available) and at z = q — n = 15 (i.e., when all the TDOAs are 
used). It is worth noticing that the availability of just a few additional TDOAs leads 
to a relevant reduction of the TDOA standard deviation, with respect to z = 0. 
This reflects also on localization, as all the algorithms monotonically improve their 
accuracy increasing the number of available measurements. Also in this case, GS 
exhibits the best accuracy after denoising, while being characterized by an unstable 
behavior using the original TDOAs. Indeed, with the original data GS improves its 
accuracy when z < 6; for higher values of z, the error bias becomes relevant and 
the overall RMSE diverges. 

It is important to highlight the practical implications of these results. Let us 
consider a microphone array composed by (n + 1) sensors. The computational 
power requested for computing all the q TDOAs is mainly due to the calculation of 
GGGs between the pairs of microphone signals. In case of limited computational 
capabilities, a typical solution would be that of computing only the nonredundant 
TDOA set. However, it turns to be more convenient to compute z additional 
TDOAs in order to fully exploit the available computational power. This enables 
better localization, without modifying the array configuration. 

7. Conclusions 

In this manuscript we reformulated the problem of source localization in the 
TDOA space. This enabled us to show that source localization has a neat interpre¬ 
tation in terms of TDOA denoising. As a simple solution to the denoising problem 
is not available, we proved that it is possible to relax the problem to a linear one, 
whose solution is based on projecting TDOAs on a linear subspace of the TDOA 
space. Moreover, we also derived the problem solution for the case in which only a 
few TDOAs measurements are available. 
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The analysis performed in this manuscript does not limit to numerical simu¬ 
lations. Indeed, each choice behind the presented algorithm is fully justified by 


means of analytical proofs that further validate and justify the works in (So et al 

As a matter of fact, in this manuscript we proved 


2008 

I and ( 

Schmidt 

1996 


that denoising has a positive effect on source localization from a theoretical per¬ 
spective. Moreover, we tested the denoising algorithm using different noise models, 
thus highlighting that the method is still valid even when noise hypotheses are not 
completely fulfilled. Finally, we also made use of different cost functions to gain an 
interesting insight on how denoising impacts on different localization algorithms. 

The extension of the relaxed denoising algorithm to the case of missing TDOAs 
have also interesting implications in real-world scenarios. As a matter of fact, it 
enables to fully exploit hardware computational capabilities in order to increase 
localization performance. As an example, by fixing the available computational 
complexity, one can tune the localization system in order to measure a given amount 
of TDOAs and fully take advantage of them in a synergistic fashion. 

According to the results of this work, it is important to note that it is possible 
to envision the development of algorithms working in the TDOA space to solve the 
complex ML source localization problem in a easier way. This will be the scope of 
possible future works. 
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